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Abstract 

We apply the orthonormalization procedure previously introduced by two of us and 
adopted in connection with coherent states to Gabor frames and other examples. 
For instance, for Gabor frames we show how to construct g{x) £ £^(]R) in such a 
way the functions gn{x) = e^"""^^^ g{x + an2) , n£l? and a some positive real number, 
are mutually orthogonal. We discuss in some details the role of the lattice naturally 
associated to the procedure in this analysis. 



I Introduction and mathematical results 



In the mathematical and physical literature many examples of complete sets of vectors in 
a given Hilbert space Ti are constructed starting from a single normalized element /o G Ti, 
acting on this vector several time with a given set of unitary operators. This is exactly 
what happens, for instance, for coherent states and for wavelets, as well as for Shannon 
systems and Gabor frames. All these examples can be considered as particular cases of a 
general procedure in which a certain set of vectors is constructed acting on a fixed element 
of H, /o, with a certain set of unitary operators, Ai, . . . ,An: fki,...,kN ■= ' ' '^aT/o, 
kj G Z for all j = 1,2, ... ,n. These vectors may or may not be orthogonal: we consider 
here the problem of orthonormalizing this set, i.e. the problem of producing a new set of 
vectors which share with the original one most of its features (we will be more precise in 
the following) and, moreover, are also orthonormal. In [1] we have constructed a general 
strategy for achieving this aim and we have applied it to coherent states, producing a new 
set of square-integrable functions which have many of the properties of coherent state and, 
moreover, are mutually orthogonal. These results will be reviewed in Section II, where 
we also describe the method. 

In Section III we apply the strategy to Shannon frame, while Section IV is dedicated 
to a deep analysis of Gabor frames. In particular we will discuss the role of what we call 
the lattice spacing in the orthonormalization procedure. 

Section V is devoted to some examples, which we will discuss in details, and to our 
conclusions. 

II Stating the problem and first results 

Let 7i be a Hilbert space, /o G a fixed element of the space and Ai, . . . , An, N given 
unitary operators: A~^ = Al-, j = 1,2, . . . , N. Let Hn be the closure of the linear span 
of the set 

^fN = {fk,,...,kr,■=A1'■■■A%^fo, h,...kNeZ} (2.1) 

which we assume to consists of linearly independent vectors in order not to be trivial. Of 
course the vectors fki,...,kN complete in Hn by definition and Hn could coincide or not 
with all of Ti. In general there is no reason why the vectors in A/tv should be mutually 
orthogonal. On the contrary, without a rather clever choice of both /q and Ai, . . . , A^, it 
is very unlikely to obtain such an o.n. set. Our aim is to discuss some general technique 
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which produces another normahzed vector if G Ti-N such that the set 



Mn = {^k„...,kr, ■■= 4'' ■ ■ ■ ku...,kNeZ} (2.2) 

is made of orthogonal vectors. Moreover, we would like this set to share as much of 
the original features of A/at as possible. For instance, if the set A/tv is a set of coherent 
states, we have shown in [1] that the new vectors ^Pki,...,kN ^^6, among the other features, 
eigenstates of a (sort of) annihilation operator, give rise to a resolution of the identity 
and they saturate the Heisenberg uncertainty relation: so they appear, in some way, as a 
set of generalized orthonormal coherent states. 

The easiest situation, = 1, goes like this: in this case we have A/i = {fk '■ = 
Jq, /c G Z} with < fk,fi >7^ ^k,i (otherwise we have already solved the problem!). 
Since A/i is complete in Tii, any element in Tii can be written in terms of the vectors of 
Afi. Let ipo G Hi be the following linear combination ipo = J^kez ^kfk, and let us define 
more vectors of Hi as = A^Lpo = J2kez fk+n = ^ fn, where we have introduced 
the operator X = Ylkez ^k^^- The coefficients Ck should be fixed by the following or- 
thogonalization requirement: < ^ni^o >= <^n,0) for n E Tj. As we have discussed in [1], 
these expansions are, up to this stage, purely formal. However, in many concrete relevant 
situations they converge and, if this is so, we find that the coefficients are 

1 e-^^ dp 

where the function a{p) is defined as a{p) = Xliez o/e*^', with aj =< A^ fo, fo >. A 
sufficient condition for a{p) to exist is that {aj} G /^(Z), which, in turns, is related to 
the nature of both A and /q. We refer to [1] and to the next Section for more results and 
examples concerning = 1. 

Let us now take N = 2. We quickly review here our previous results on coherent 
states. 

Let q and p be the position and momentum operators on the Hilbert space H = £^(]R), 
[q,p] = iH, and let us now introduce the following unitary operators: 

U{n) = e''^('^i9-'^2P), D{n) = e'^''^-^^\ Ti := e^"^ Ta := e'^^P. (2.4) 

Here a is a real constant satisfying = 27tL for some L G N, while Zn and b are related 
to n = (ni, na) and q, p via the following equalities: 

z„ = -^(na + mi), b= ^{q + ip). (2.5) 
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With these definitions it is clear that 



U{n) = D{n) = (_l)i"i«2Ti"^T2"^ = (^_l^Lmn2j.n2j.ni^ ^2.6) 

where we have also used the commutation rule [Ti,T2] = 0. 

Let ipo be the vacuum of b, b(po = 0, and let us define the following coherent states: 

(/.(f) := Tf^T^^Vo = T^'T^^^o = (-l)^"^"^t/(n)<^o = (-l)^"^"^Z^(n)<^o, (2.7) 

[21 [3]. It is very well known that the set of these vectors, C*-^^ = {^pi^\ n G Z^}, satisfies a 
number of relevant properties l]. However, it is also well known that they are not mutually 
orthogonal. Indeed we have J^^'' :=< Lp^n\^Q >= (— l)^'^i"2 g-f which is only 

nearly zero if L is large enough and (rii, 77-2) 7^ (0, 0). 

Our aim is to construct a family of vectors S^^'^ which shares with C*-^-* most of the 
above features and which, moreover, is made of orthonormal vectors. We have shown in 
[1] that this is possible, in suitable Hilbert spaces, if L > 1, while complications arise for 
L = 1. Incidentally, we have shown that the set C^^^ is complete in Ti, if and only if L = 1. 

Let us now define, for each L > 1, the following set: 

/ii := linear spanj^J^f 77, G Z^l . (2.8) 

We know that hi = Ti, while, for L > 1, hi G Ti. It is further clear that hi is an Hilbert 
space for each L, since it is a closed subspace of 7i. 
Now we define 

:=E^'4''b (2-9) 

fc6Z2 

Of course this means that \l/g^^ := Xlfeez^ '^i^'' V^i^'* ^^^^ because of the commutativity of 
Ti and T2, we also have 

^(L) = J^n, j.n2^{L) _ ^2.10) 

Therefore the new set constructed in this way, := {^^^\ n G Z^}, is invariant 

under the action of Ti and T2, exactly as the set C^^\ independently of the choice of the 
coefficients of the expansion c^^"*. These coefficients, however, must not be chosen freely: 



"'^(a) C^^-' is invariant under the action of TJ*^, for all nj, j = 1,2; (b) each ipi^"^ is an eigenstate 
of b: htp\t^ = Zn'J'n^', (c) they satisfy the resolution of the identity on a certain Hilbert space /il, see 
Stigz^ Wn^^ X Vif^l = li: where Ul is the identity on /i^; (d) They saturate the Heisenberg 
uncertainty principle: if we call (AX)^ =< > — < X >^ for X ~ q,p, then Aq Ap = i. 
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they are fixed requiring that the vectors in the set £^^^ are orthonormah < \I/n , >= 
6n,s- This will fix (not uniquely!) the value of the ci^'''s, with a procedure which extends 
formula fl2.3l) . These coefficients, if they exist, can be found as follows: 

cf^ = ^ / / ^ " dP, (2.11) 

where Fl{P) = Emez^ (-l)^'"^"'^ g-f L(mf+"^i) e'^'^. 

As already mentioned, in [1] we have seen that L = 1 and L > 1 are really different 
situations: if L > 1 the above formulas are well defined and produce a set £^^^ which has 
the same properties of the vectors in C^^^ and, moreover, is made of orthonormal vectors. 
On the contrary, if L = 1, the procedure does not work: the series are not convergent and 
Fl{P_) has a zero in [0, 27r[x [0, 27r[. This is not related to our method but it is rather a 
well known feature of coherent states: indeed it is a standard result in functional analysis 
that in Ti. no orthonormal set of coherent states can be constructed, [2]! 



Ill Frames of translated 

We begin this section by recalling few known results concerning frames of translated. 
After that, we will use our method to construct an o.n. set of translated. 

Assume that G £^(]R) and consider the set of the form JF = {(/)(■ — k)}kez = {Tk4>} 
where is related to the translation operator in 02.41) as follows: := with a = 1. 
Christensen, Deng, and Heil in [1] by using the Beurling densities proved that JF cannot 
be complete in all of £^(M) and that it can, at most, be a frame for a proper subspace 
of /^^(M). This is in line with what we have discussed in the previous section: indeed the 
closure of the set A/i was not required to be all of H. What our strategy can produce here 
is, at most, an o.n. set spanning the same Hilbert space which is spanned by the set JF. 

Let us consider the function 

$:M^M, $(p) = ^|0(p + A;)|2 (3.1) 

kez 

where cj) is the Fourier transform of 0, 0(p) = 't'i^) e~*^^ dx. The function $ is 

1-periodic and belongs to >C^([0, 1]). This function can be used to check wether the set 
is or not a frame of translates. Indeed, see for instance [5] and references therein, the 
following proposition holds true: 

Proposition 1 Let (p G £^(]R) then for A, B > the set JF is an [A, B)-frame if and 
only if 

A < ^{p) < B 
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a.e. inp e [0, 1] \ C, where C = {p e [0, 1] : = 0}. 



Remark: we want to stress that we are here adopting a notation which is shghtly 
different when compared with [5]: in particular we still call frame a certain set even if 
it only spans a subset of the original Hilbert space, at least when this aspect is clear 
and does not lead to confusion. In other words, whenever the situation is clear from the 
context, we will say that a certain set of translated is a frame even if it does not spans all 

of n. 

Other results on frames of translates are discussed in many details in [5], to which we 
refer for more details. 

Let us now apply our procedure for = 1 to a set of translated. For that we suppose 
that an (A,B)-frame Afi = {(pj G 7i}jez is generated by a single fixed vector ip and a 
single unitary operator T = T2: 

ipj = T^if, j G Z. 

So there exist A,B > such that 

where h is the subspace of H spanned by the vectors of Af. The construction of the o.n. 
set goes as in the first part of the previous section. Therefore we can define the coefficients 
ttj as usual, aj =< T^ip,ip >, and from these we find the function a{p) and finally the 
coefficients Cfc's which produce the vector ipo = Ylik& ^kfk as in (12.31) . 

An example of this construction was already essentially discussed in [1]. We consider 
this same example here from a slightly different perspective. 

Example:— Let g{x) = X[o,a[i^) characteristic function in the interval [0,a[, 

with 1 < a < 2. Of course we have 

Ml = {gn{x) := T^'gix) = X[n,n+a[{x), U E Z} . 

The overlap coefficients aj can be written as aj = aSj^ + (a — 1) {Sj-i + so that 
a{p) = a + 2(a — 1) cos(p). This is a nonnegative, real and 2-7r-periodic function and is 
never zero in [0, 27r[ since it has a minimum in p = vr and a;(7r) = 2 — a > 0. If we fix, as 
an example, a = |, we can compute analytically ^^g^ |qP = J^^ = We refer 
to [T] for further considerations on this example. Here we just want to comment that the 
function = J2kez IdiP + (13. ip can be written as = a + 2(a — 1) cos(27rp). 
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and this satisfies the inequahty 2 — a < ^{p) < 3a — 2. Therefore, using Proposition 1 
above, we conclude that the set {gn{x), G Z} is a (2 — a, 3a — 2)-frame (which, of course, 
spans a proper subspace of £^(]R)). So we have a frame which, by means of our strategy, 
produces an o.n. set which spans the same Hilbert subspace of £^(M) and which is stable 
under translations. 

We would also like to mention that a different orthonormalization procedure is well 
known to people working on wavelets since it may be useful in the very first steps of a 
multi-resolution analysis to construct an ortho normal basis of a certain subspace of 
Vq, starting from a given Riesz basis of Vq. This technique, which is reviewed in [6], is 
mainly based on the following fact: in (f{x), together with its integer translated (f{x + n), 
n G Z, is a Riesz basis in Vq C £^(M), then the inverse Fourier transform (p{x) of the 
function (Xliez + 27r/)|^) ^ (^(p) is such that the set {<^(x + n), n G Z} is an 
orthonormal basis of the same set Vq. 

It is clear that this assumption exactly coincides with our original hypothesis: any 
Riesz basis is a set of vectors {Tk4>}kez which are linearly independent and a system of 
generator in the Hilbert space they generate, by definition. 

It is also worth mentioning that the two techniques both rely on Fourier analysis and, 
from this point of view, they look similar. Nevertheless they are different and inequivalent 
since they produce, in general, different results starting from the same seed function g{x) 
and, moreover, since the difficulties of the computations are not necessarily comparable. 
For instance, if we consider the previous example, our orthonormalization procedure pro- 
duces, at the end, the following function: ( tt- fn j — - '^^ I Yin n+alix), and 

' ' y2n JO ^a+2{a-l)cos(p) J ^1"'"+"!^ ^' 

the integrals can be easily computed (Remember that we need to compute only very few 
of them, because general arguments on Fourier series show that the coefficients c„ go to 
zero very fast with n, see below). A bit more difficult, in our opinion, is to compute the 
function if{x) using formula above. However both these approaches in general can only 
be applied within an approximation scheme. Finding an exact solution is, but for few 
examples, quite difficult if not impossible. From this point of view our strategy seems 
convenient with respect to the other one since it is very natural to construct a perturbative 
approach. Indeed we only need to check if the function a{p) has some zero in [0, 27r[. If 
this is not so, in fact, we know that the coefficients Cm goes to zero faster than any inverse 
power in |m| and, therefore, the series defining the new o.n. set converges very fast. This 
means that only the first few c^'s are required to get a reasonable approximation of the 
solution we are looking for, and this can be easily done numerically. 

Remark: It is maybe not surprising to notice that both in our and in the above 
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approaches the following integrals play an important role 



ak= 4>{t)(t){t + k)dt. 



J —oo 



This is because is a measure of the non orthogonality between the original functions: 
the smaller is the difference of from 6kfi, the closer is the set + /c), G Z} to an 
orthonormal set. 

IV Gabor frames 

IV. 1 What we get from (/c, g')-representation 

We begin this section by recalling few results on the generalized (/c, g)— representation 
which we have introduced in [1], and which extends analogous results discussed, for in- 
stance, in [7]. 

Let T(a) = e^^"" and r(6) = e*^'', with ab = 27rL, for some natural L. It is clear that 
for all possible L G N the two operators still commute: [T(a),r(6)] = 0. For each given 
A > let us define the set of (generalized) functions 




(4.1) 



where (k, q) G := [O, ^ [ x [0, a[. If is the generalized eigenvector of the position 
operator g, q^x = x^x, we write ^i^g''\x) as ^i^g''\x) =< ^[^f^ >. 



Proposition 2 With the above definitions the following statements hold true: 



1. 



(4.2) 



2. 




(4.3) 



3. 




(4.4) 



where the usual Dirac bra-ket notation has been adopted; 



[ ^tf{x)^[^:;){x)dx = 5{k-k')5{q 




(4.5) 
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The proof of these statements does not differ significantly from the standard one, and 
will be omitted here. We just want to remark that, for general a and a', we find that 
T(a)$[.'^'" 7^ e^^^^^^^°'\x). In other words, in general <^^^^'^\x) is not an eigenstate of 
T(a) if a 7^ a' . 

Also, it should be noticed that the value of the parameter h entering in the definition 
of t(6), is fixed by requiring that T and r commute but play no role in the definition of 
the lattice cell which on the other way is defined by an extra positive parameter, A, 
which needs not to be related to h itself. However, quite often in applications A coincides 
with a and with h. 

We now use this generalized (fc, g) -represent at ion in connection with Gabor frames and 
in the attempt of getting an o.n. set of functions in £^(]R). The starting point is not very 
different from that in Section II: let q and p be the position and momentum operators on 
the Hilbert space U = C^iR), [q,p] = zl, and let U{n) = e^"("i«-"2P), D{n) = e'^^^-^^\ 
Ti := e*"''^ and T2 := e~*"^ as in (12.41) As before, a is a real constant such that = 27rL 
for some L G N, while Zn and b are related to n = (rii, 71.2) and q, p as in (12.51) . As in (12. 6p 
we have U{n) = D{n) = (_i)i^«i"2 j^j^i j^^"2 ^ ^_^^Ln,in22.n2 ^^^^^ ^^^^ ^^^.^ ^^^^ 

the commutativity [Ti,T2] = 0. 

Let now g{x) be a fixed function in Ti, and let us define the following functions: 

gl^\x) := T^^T^- g{x) = T^-T^^ g{x) = (-l)^"^"^f/(n) ^^(x) = 

(4.6) 

We call C(^) = {gn (x), n G Z^} the set of these vectors, which is the standard set of 
coherent states if g{x) is just the vacuum of the operator b. These functions, which quite 
often in the literature are written as e^'^^'^°^ g{x — n2XQ) for some positive Xq and ujq, are not 
mutually orthogonal for a generic g{x). For instance, they are not orthogonal for coherent 
states. Here we will assume always that the functions in C-^-* are linearly independent. 

Our aim is to construct a family of vectors £^^^ which shares with C*-^-* most of its 
features and which, moreover, is made of orthonormal vectors: we are specializing our 
original settings of Section II to = 2. Of course, we are interested here in considering 
a starting function g{x) such that < gl%'\g >^ 5n,o- Also, we would like to know if for 
some choice of g{x) G £^(]R) and L > 1 it is possible that C^^^ is a complete set in £^(M). 
In other words, if we define Tig.L := linear span {gii'\x), nE Z^}, we wonder whether a 
clever choice of g and L produces Ti-gx = £^(M). These questions can be answered rather 
easily by making use of the generalized {k, g) -represent at ion. 

We first remark that, introducing the generalized common eigenstates of Ti and T2 as 
in the beginning of this section, we have Ti^^f^""^ = e*'J''$^^^'"-' and T2$[j^g'"'* = e^^"^^\f^"'\ 
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Therefore < gnK^i'^f^ >= e"*"!'^'' g-^^^a^ 
tion of the identity in (I4.4I) . imphes that 

>'2t7/A 



< g, ^^"^'"^ > which, in turns, using the resolu- 



9,m 



--< g[!'\g >-- 



dke 



-ikn2A 



(4.7) 



where we have introduced the {k,q) representation of the function g{x) as g{k,q) 



--< 



'^fei/'^^fi' ^- ^^^^ ^'^^ that, as in [T], the role of L is crucial here. 



k,q 



If L = 1 the set S'^^'^ = {e-^'^"^^ e"^""!'', (^1,712) G Z} is complete in ^^^D^^)), so 
that the functions in C^^^ can be mutually orthogonal if and only if g{k, q) is a phase so 
that \g{k,q)\ is constant almost everywhere (a.e.) in This is the reason why the 

coherent states are not mutually orthogonal: the (/c, q) representation of the vacuum of 
the annihilation operator h is not simply a phase! 

If L > 1 the conclusion changes. Let us consider, as an example, L = 2. In this case the 
set iS*^^-* is no longer complete in £^(n^^)) so that equation (14. 7p does not automatically im- 
plies that \g{k, q) \ must be constant. On the contrary, splitting the double integral on O^^'^ 
in two contributions over {k, q) G [0, 27r/yl[x [0, a/2[=: □[^g and (/c, q) G [0, 271 / A[x [a/2, a[, 
we can write as 



/(2) 



2-k/A 



dk e 



-ikn2A 



a/2 



rfge-^^"^'^ {\gik,q)f + \gik,q + a/2)\') 



which is equal to (5„^o if and only if \g{k, q)\ +\g{k,q + a/2)\ is constant almost everywhere 
yl since 5^^^ is complete in £^(n^^2 



for {k,q) G □^'^2 since 5^^^ is complete in C'^{D^^^). This result is easily generalized: 



Proposition 3 The setC^^^ generated by a given square-integrable function g(x) is made 
of orthonormal functions if and only if in the {k, q) -representation g satisfies the following 
equality almost everywhere 



\g{k,q)f+ g {k,q + 



L 



+ 



g \k,q + 



2a 



+ ■ 



g {k,q 



(L- l)a 



LA 
2Tia 



(4.8) 



Of course, once such a g{k, q) has been found, then the related g{x) can be simply obtained 



as 



g{x)= [ [ dkdq^[^f{x)g{k,q) 



(4.9) 



Any g{x) which cannot be written in this form cannot produce an o.n. set C^^\ As a 
consequence, this is exactly the kind of functions we will consider in the rest of this section, 
since we are interested here in producing an orthonormal set starting from a different set 
C^^) which is not orthogonal from the very beginning. 
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The generalized {k, q) representation can also be used to discuss the problem of the 
completeness of C*^^^ in £^(]R). Before doing this, however, it may be worth observing that 
by definition C^^^ is complete in 7ic,,L but, since in general 7ic,,L is only a subset of £^(M), 
we have no information about the completeness of C^^-* on this larger Hilbert space. In 
other words, if we find conditions for C^^-* to be complete in £^(R) we also find conditions 
for 'Hg^L to coincide with £^(]R). 

To answer these questions we first recall that the set C^^^ is complete in £^(M) if, 
given a square integrable function h{x) which is orthogonal to g^\x) for all n^l?'^ then 
h{x) = a.e. in M. Using the properties of the functions we can write 

JO Jo 

where, as usual, g{k,q) =< ^k^g"'\g > and h{k,q) =< ^\fg"'\h >. Again, to find 
conditions for these scalar products to be zero Wn E 1? , it is convenient to consider 
separately the two cases: L = 1 and L > 1. 

If L = 1, using as before the completeness of the set S^^^ in /^^(D^"^^), we deduce that 
< gil'\ h >= y n G Zi^ if and only if g{k, q) h{k, q) = a.e. in Therefore, if g{k, q) 

is zero at most on a set of zero measure, this implies that h{k^ q) =^ a-e. in □'^"^^ and, as 
a consequence, that h{x) = a.e. in M. 

Remark: the vacuum of b, which generates the set of coherent states, satisfies this con- 
dition, [8], and therefore the related set C*-^-* is complete in £^(M). Examples of functions 
which do not satisfy this condition can be easily constructed. Consider, for instance, a 
vector g which in the (/c, g)-representation is equal to 1 for {k,q) G and zero other- 
wise. In the ^-representation this function looks like g{x) = X[o,a/2[{.x), where 
is, as usual, the characteristic function of the set /. For what we have just shown this 
function does not generate a complete set in £^(M) since non zero functions h{k,q) for 
which g{k, q) h{k, q) = a.e. in O^^^ can be easily found. 

Let us now consider what happens for L > 1 and, to be concrete, let us fix L = 2. 
Splitting the integral as we have done before, we deduce that < gil'\ h >= for all nEl? 
if and only if is zero the following combination g{k, q) h{k, q) + g (^k, q + ^) h {k, 9 + f ) • 
In other words, < g'^\ h >= for all n G if and only if 

g{k, q) h{k, q)+g(^k,q + ^^h[k,q+^^ =0 a.e. in dJ^^ (4.10) 

and, as a consequence, C^^^ is complete in for some given g{k,q) if and only if 

the only solution of equation fl4.10p is h{k,q) = a.e. in n^g- course, any function 
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g{k, q) which is zero on a set P C Di/^ of non zero measure cannot produce a complete 
set, because equation (14.101) would have a non trivial solution: it is enough to choose a 
function h{k, q) which is zero only outside P! Hence, in order to obtain a complete set, 
it is surely necessary to assume that g{k,q) is zero at most on a set of zero measure. 
However, even under this assumption it is easy to check that other non trivial solutions 
of this equation do exist, and therefore C^^^ cannot be complete in £^(M). Let indeed 
So{x) be a fixed function in £^(M) and let fJ,{k,q) := e*^"/^ Zlnez e*"'^^So(g + na/2). 
fi{k, q) belongs to C'^{0^yl) since J J^(a) \fi{k,q)\'^ dkdq = ^ |so(x)p (ix, which is finite 

1/2 

because So(x) G £^(M). Using now the boundary conditions g{k, q + a) = e*^"^ g{k, q) and 
fi{k,q + a/2) = — e*^^ g), it is easy to check that h{k,q) := g{k,q + a/2) fi{k,q) 
satisfies equation fl4.10p and is different from zero a.e. if So{x) is chosen, e.g., positive, 
since in this case g{k, q) a.e. in O^yl- 

Let us now summarize the above results: the set C^^^ consists of o.n. functions if and 
only if g{k,q) satisfies equation (^7^. Moreover, if g{k,q) ^ a.e. in D^^^, then C^^^ is 
complete in £^(M). Finally, if L > 1, there is no choice of g{k, q) for which the set C'-'^^ is 
complete in £^(]R). Therefore, for any given g and for all L >2 we have 'Hg^L C £^(]R). 

This result reflects, in a certain sense, what happens for frames of translates which 
are surely not complete in all of £^(M), as we have discussed in the previous section. 



IV. 2 Back to the orthonormalization problem 

We are now ready to apply our orthonormalization procedure to the set gn{x) in (14. 6p . 
with the only requirement that C'^^'^ must consists of functions which are not orthogonal 
from the very beginning and linearly independent. As in [T] and in Section II we define 

4'n^):=E (4-11) 



which implies that '^^q\x) := Xlfeez^ '^i^'* 9l^\^) ^^d, because of the commutativity of Ti 
and T2, that 

^^^Hx)=T^^T^-¥l\x). (4.12) 

Therefore the new set constructed in this way, (x), n G Z^}, is invariant 

under the action of Ti and T2, exactly as the set C^^\ independently of the choice of the 
coefficients of the expansion c\^\ These coefficients ci^'' are, as usual, fixed (non uniquely) 
by requiring that the vectors in the set S^^^ are orthonormal, < \E'^^\ \E'i^'' >= 6n,s, or, 
equivalently, that < ^'^^ \E'o^^ >= 6n,o. Indeed, if they exist, the coefficients can be found 
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27r /.27r ^-iP-k 



as before: 

where Fl{P_) = ^^^12 Im^ e*— — . As in [T], the behavior of the coefficients in (14. 131) is 
directly related to the nature of the convergence and to the zeros of Fl{P_), which, in turns, 
are related to the coefficients of the overlap between gm\x) and g{x), Im^ =< gm\g >■ 
In particular the following results immediately follows from basic facts in the Fourier 
series analysis and from the results in [1]: 

1. if J^f^ e then Fl{P) is continuous, non negative, {2n, 27r)-periodic and 
bounded; 

2. if e s(Z2) then Fl{P) is a non negative function; 



3. if /^^ are such that E^gz^ < 1 then Fl{P) ^ for all P G [0, 27r[x [0, 27r[ 
and, as a consequence, {cm''} G s(Z^); 

4. if G s(Z^) then the sequence \I/^^]v '■— J2k-\\k\\<N ^^k^'d^k^ converges in the ||.||2- 
norm so that its limit, '^^n \ belongs to 'Hg^L- 

If we introduce the operator Xl = Emez^ ^ii^ '^P then ¥n^ = Xl g^n \ It is also 
possible to extend the following proposition, originally given in [1]: 

Proposition 4 Suppose that {c^^} G /^(Z^) and that Fl{P) ^ OforallPe [0, 27r[x [0, 27r[ 
Then is complete in 'Hg^L if o.nd only if X^ admits a bounded inverse. 

The proof is very close to that given in [1] and will not be repeated here. 

It is rather interesting to observe that the same sum rule which was proved, under the 
assumptions of this Proposition, in [1], that is E/cgz^ '^i^^ '^i^^ — 1; ^^^o holds true in the 
present settings. Here the a\^^ 's are the coefficients of the expansion of the non-orthogonal 
functions g^^ in terms of the orthogonal ones, g^\x) = EmGZ2 (^ir '^l5!-m{^)- 

Explicitly we find ai^^ = J^"" q-^p-ir ^ Fl{P}. Again, we don't give here the 

proof of this result which is based on the Poisson summation formula since it is very close 
to that given in [1]. We just want to stress that for each starting function g{x) we get a 
different sum rule. 

The vectors of the set 8^^^ have the following properties: 

E 14'^ X = E X = h^L, (4.14) 



riL) 
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where Ig^L is the identity operator in Tig^L, which, for what we have proven at the beginning 
of this section, can be the identity in £^(M) only if L = 1 and g is taken conveniently. 



same eigenvalue. 

Here as in [1] the role of L is very important and in particular we are not sure a priori 
that Xj^^ exists for all values of L. For instance, if L = 1 it is possible to adapt the 
same reductio ad ahsurdum argument given in [1] to conclude that our orthogonalization 
procedure cannot work. More in detail, suppose that we have constructed the o.n. set 
£W which is a basis of £^(M). Then we conclude, as in [T], that the original set C^^^ is an 
o.n. basis of £^(M) as well, which is false. 

Before considering some examples a final remark is in order: for Gabor frames it 
is widely discussed in the literature, see [B] for instance and references therein, that 3 
different regions appear in their analysis and they give rise to different mathematics. 
More in details: the set {gnfx) = e^°'^^^g{x + an2), ni,n2 G Z} can be an orthonormal 
basis only if = 27r (but with poor localization properties). It can be a frame for all 
of if < 27r (and it may have good localization properties). Finally, it cannot be a 
frame for all of ?i if > 27r. Since our strategy works for = 27rL, L = 1,2,3, .. ., 
this automatically excludes the case < 27r. This constraint, in our approach, is due to 
the important requirement that [Ti,T2] = which we hope to weaken in a close future. 
However, as we have widely discussed here and in [1], already for > 27r we find a very 
rich mathematical setting, which produces many interesting results. 

V Examples and conclusions 

The first example of this construction produces the classical coherent states. We don't 
discuss this example here since few results have already been given in Section II and a 
complete analysis can be found in [1]. In this paper we concentrate our attention on other 
examples whose computations can be carried out almost completely analytically. Other 
examples can be discussed in a totally analogous way. 
Example 1. 

Let us consider the following normalized compactly supported function. 



Also, if by chance gn (x) 
coherent states), then '^if^\x) 



is an eigenstate of some operator ki (as it happens for 
is an eigenstate of the operator Kl = Xl Xf-^, with the 
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It is clear that the different gn{x) = e^"""^^^ g{x + an2) are not all automatically orthogonal, 
even if < gn^go > is surely zero if \n2\ > 2 because, in this case, the supports of the two 
functions have empty intersection. It is easy to check that, choosing L = 4 so that = Stt 
to simplify the computations, Fl{P) = 1 + | cos(P2) =: Fl{P2) which is independent of 



Pi. As a consequence, the coefficients c\^^ in fl4.13p can be written as c)^^ = ^fci.oCfca, 



where 



1 

2^ 



27r g— iP2fc2 



dPo. 



'0 ^/F;XP2) 

These coefficients can be easily computed: Cq = 1.11308, c±i = —0.216769, c±2 = 
0.0625106 and so on. Needless to say, because of the analytic properties of -Fl(P2) which 
is never zero in [9, 27r[, they tend to zero faster than any inverse power of k2 so that all 
the series which appear along our computations are surely convergent (in the strongest 
topology). 

From formula (14.111) we find that \E'(x) = X]fc^fcfi'(^ + ^^)- P^^^ following 
figure \E'Ar(x) = ^^_7v Ckg{x + ka) for different values of A^. 



Figure 1: ^'(x) for iV = 0, = 1, = 2 (first row) and A^ = 3, A^ = 4 (second row) 



It is just a technicality to check that the overlap between '^{x) and T"^T^^\E'(a;) can 
be written in terms of the coefficients as follows: 

< ^(f), ^ >= [q(q+„, + (q+„,+i + Q+„,_i)/3)] 5n,,0 (5.1) 

and we find that, replacing Ylii^i with Yld=-N (15.11) . for A^ = 3 we have H^&lP = 1.00001 
while the modulus of the overlap between two different wave-functions does not exceed 
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0.00605858. Therefore we conclude that this is aheady a very good approximation, which 
can however be improved if we take = 4: in this case we find = 1 and the modulus 
of the overlap between two different wave- functions is always less than 0.00208293. 

The computation of directly follows from its definition: Xi = ^„2ez^"^^ -^2"^- 
inverse, X£^, can be computed in complete analogy: since g^\x) = Xlmez^ '^ir ^n+mi^) = 
X-'¥r!:\x), we get Xj:' = Y.n,&^n^ where = ^ ^ dP e-'^"- ^/F^) . We 
find, for instance, = 0.96857, = 0.175095, d^^s"^^ = -0.016400, and so on. 

Example 2. 

Let us consider now the following function: 



9[x) 



Nb exp (^azp) , \x\<b 
0, otherwise. 



where b = a = 2'kL, has been introduced only to simplify the notation. A*";, is a 
normalization constant which depends obviously on b and, as a consequence, on L. Once 
again we consider a function with compact support just to simplify the computations, since 
with our choice < gn,g > is zero if \n2\ > 2. Hence, the function Fl{P_) can be written 
as Fl{P) = /o(Pi) + 2A(P0 cos(p2), where /o(Pi) = En.ez ^n^^''^'''' ^nd j\{P^) = 
SniGZ -^nt^i • '^^^ simplest results are obtained when both these functions can be 
reasonably well replaced by their constant main contributions, that is those contributions 
coming from rii = 0. Indeed, if this can be done, then we get Fl{P_) ~ 1+2 |) cos(P2) =: 
Fl{P2) which does not depend on Pi, so that, once again, c\^^ = ^A.-i.oCfc^'', where c[.^^ = 
^ In^ ~lr==dP2 and "^^^^x) = ^uC^'^ g^^\x + ka). Of course, replacing /o(Pi) and 

/i(Pi) respectively with their first contributions I^^J = 1 and Iq^^ is possible only if we 
have some control on what we are neglecting. However, if L = 1, it is easy to check 
that -^(d^l) is smaller than that part of /2(P2) we are neglecting, 5/2(P2), so that replacing 
Fl{P_) with Fl{P2) is a dangerous operation! This will be reflected on the fact that, as 
we will see, the orthonormalization procedure does not work in this case: we get a new 
family of functions which are not all mutually orthogonal as we would like. On the other 
hand, if we consider L = 2, then |5/2(P2)| < l-^o^"*! but they are of the same order of 
magnitude: we have some chance that the procedure works but we cannot be sure at this 
stage. For L > 3, finally, |5/2(P2)| ^ l-^o^^l ^^^^ we expect to get a good orthonormal 
set for Tig^L- 

These are exactly the conclusions that we obtain at the end of the computations: if 
L = 1 the function ^^^\x) = J2k^k^^ d'^^^i^ + ^'^) ^^^^ ^^^^ — 0-96, while, for 

instance, < ^f*^^), >~ 0.1. We see that the related functions do not produce an 
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orthogonal set. 

Let us see what happens if L = 2. In this case we find II ^(2) II 2 ^ 0.99996 and 
< ^(^), >^ 0.0001, which are quite promising. However, we also find that < 

v|/{2)^ v[/|^''^^ >~ 0.037 which shows that for L = 2 the approximated \l/^^^'s constructed 
as shown above are orthogonal up to corrections which are small but not too much. 
However, as expected, this out of orthogonality parameter gets smaller and smaller when 
L increases: already for L = 3 we find < ^(i o) —0.012, and all the other scalar 
products are much smaller. The norm of "^^^^ is essentially 1. In the following figure we 
plots the different approximations of "i/^^^x) for L = 3. 




Figure 2: ^(x) for iV = 0, iV = 1, = 2 



Example 3. 

Let us consider now the following function: 



^ I COS (t) , \x\<a 
I 0, otherwise. 



Also in this example we consider a function with compact support since it allows to 
perform almost all the computations analytically. As before < gn,g > is zero if |n2| > 2, 
and this makes the computation of Fl{P_) simple. Indeed, if we introduce the function 
f{Pi) = Em.ez then we get Fl{P) = 1 + 2 f (P,) cosiP^) . Of course, due to 

this formula, c^f'' is not the product of 5„j^o and a coefficient which only depends on 
n2. Nevertheless, this is exactly what happens if we consider the first approximation of 
/(Pi) ~ ^, in analogy with what we have done in the previous example. Again, this 
replacement is justified especially for values of L larger than 2 and we find c^^"* = 6ki,o c^j, 
where c,, = ^ dP^, F^iP^) ^ 1 + f cos(P2), and ^{x) = ~g{x + ka). 

The coefficients can again be easily computed: cq = 1.0997, c±i = —0.20105, c±2 = 
0.0545131 and so on. Of course, because of the analytic properties of Fl{P2), they tend to 
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zero faster than any inverse power of k2 so that even in this example all the series which 
appear along our computations are surely convergent. 

We plot in the following figure \1/Ar(a:) = J2k=-N '^kOix + ka) for different values of 
N and for L = 1. Not many differences appear in the shapes of the functions for larger 
values of L. The only major difference is a bigger support. 




Figure 3: ^(x) for iV = 0, iV = 1, iV = 2 (first row) and iV = 3, iV = 4 (second row) 



For what concerns the check of the orthonormality, we first notice that formula (15. ip 
does not hold here and should be replaced by the following one, which can be easily 
analitically derived: 



< 



l&Z 



+ 



1 



7r(l - 4nfL2) 



^ Q (q_„2+i + 



(5.2) 



As we see, this formula explicitly depends on L, as it should. Also, replacing the sum 
on / G Z with a finite sum on / = — A^, . . . ,N, we can compute the norm of "^{x), which 
already for A^ = 4 is equal to 1 but for an error smaller than O(10~^). An interesting 
result concerns the overlap between two different functions and \&(x). If rii = we 



find that 



< ^1^^ ^ > 

^ ^0,n2' ^ 



< 10 already for A^ = 4, while if A^ = 5 it is even smaller 



than 5 X 10 ^, V?7,2 G Z, independently of L. However, L plays a role in the computation 



of 



(L=l) 



if rii 7^ 0. Indeed in this case we can easily check that, for instance, 

.155. This result shows that, at it was already discussed, our procedure 
cannot produce an on basis when L = 1. A different conclusion is obtained if L > 1. 
Indeed, already for L = 2, we find that the maximum overlap is given by < ^* > 
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.031, already for = 4. This is a reasonable approximation which gets better and better 
as L increases: the same quantity equals 0.013 if L = 3, 0.007 if L = 4 and so on. To 
improve the approximation we should consider a better approximation for /(-Pi). However, 
such an improvement necessarily produce a more complicated expression for c^'' , and will 
not be considered here. 

More examples are easily constructed starting from functions with compact support 
in [—a, a]. For these functions the analytic expressions of Fl{P_) can be found with no 
particular difficulties, and reasonable approximations can also be deduced in most of the 
cases. The situations becomes technically more difficult when the starting function g{x) 
has no compact support, as in [T]. In this case all the computations are usually more 
delicate even if, in principle, they still produce an o.n. basis in the Hilbert space Hg^i- 
More than constructing other examples we are interested in extending the orthonormal- 
ization procedure to a slightly different situation, i.e. to the case in which the original 
set of functions {T"^ T^^5f(x), ni,n2 G Z} is constructed using two unitary operators Ti 
and T2 which do not commute. This analysis, in fact, would produce interesting outputs 
related to Gabor systems for all of £^(M) and to systems of wavelets. 
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